Three-dimensional fast electron transport for ignition-scale inertial fusion capsules 
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Three-dimensional (3D) hybrid PIC simulations are presented to study electron energy transport 
and deposition in a full-scale fast ignition configuration. Multi-prong core heating close to ignition 
is found when a few GA, few PW beam is injected. Resistive beam filamentation in the corona seeds 
the 3D current pattern that penetrates the core. Ohmic heating is important in the low-density 
corona, while classical Coulomb deposition heats the core. Here highest energy densities (few Tbar 
at 10 keV) are observed at densities above 200 g/cm 3 . Energy coupling to the core ranges from 20 
to 30%; it is enhanced by beam collimation and decreases when raising the beam particle energy 
from 1.5 to 5.5 MeV. 



PACS numbers: 52.57.Kk, 52.65.Ww 

Fast ignition of inertial fusion targets [H, 0] offers a 
promising alternative to the standard scheme of central 
hot spot ignition 0, Separating fuel compression 
from hot spot heating is expected to reduce compression 
and symmetry requirements significantly. Here we study 
Tabak's original proposal |l| to use a laser-driven elec- 
tron beam to ignite the compressed core. Atzeni 0, [f| 
estimates that an energy of about 20 kJ is required to 
ignite 300 g/cm 3 Deuterium- Tritium (DT) fuel when de- 
posited in 20 ps on a 20 /im spot radius with 0.6 g/cm 2 
stopping range. This implies a giga-ampere (GA), peta- 
watt (PW) pulse of 1 MeV electrons. Little cones may be 
used to generate such beams and to guide them through 
the plasma corona 0, Q. Recent experiments on cone 
guiding with 300 J, 0.5 PW laser pulses have demon- 
strated excellent (20 -30%) energy coupling from laser to 
core [6, 7]. For full-scale fast ignition, a 100 kJ, multi- 
PW beam is required including energy coupling. It will 
carry a few GA current, 10 times larger than the Alfven 
current that may limit transport due to magnetic self- 
interaction. Here the beam has to be transported over 
a distance of 100 - 200 /im between cone tip and core 
through a high-gradient plasma profile. This is the topic 
of the present paper. 

In plasma, the beam current is compensated by return 
currents, which suppress the magnetic fields. But this 
beam is subject to filamentation instability. For colli- 
sionless plasma, linear growth rates have been studied 
in 0, [l(| , and particle-in-cell (PIC) simulation was used 
to trace the nonlinear evolution [ll|, [12]. This applies 
to low plasma densities, comparable to the beam den- 
sity (ribeam/npiasma > 0.1) but not to the high plasma 
densities considered in this paper. Full-scale PIC sim- 
ulations are not yet feasible. They should include col- 
lisions and plasma resistivity to properly describe the 
return currents. Here we use a hybrid model adequate 
for describing self-magnetized transport in high-density 



fuel. It treats only the relativistic beam electrons by PIC 
and models the background plasma by the return current 
density j r , tied to the electric field E—rjj r by Ohm's law 
with resistivity rj. Maxwell's equations are used in the 
form VxB=^oj and VxE=— dH/dt, where j=jh+j r is 
the sum of beam and return current density. The dis- 
placement current and charge separation effects can be 
neglected since in this high-density environment relax- 
ation times and Debye lengths are much smaller than 
the sub-picosecond and micrometer scales of the resistive 
filamentation investigated here. The beam deposits en- 
ergy into plasma electrons in two ways: by direct classical 
Coulomb deposition and via return current ohmic heat- 
ing with power density ?yj 2 . Electrons and ions are cou- 
pled by thermal energy transfer. A plasma density pro- 
file constant in time with equal electron and ion number 
densities is assumed. This model was pro posed by Bell 
[1 31 ] and further developed by Davies [lj] and Gremil- 
let et al. [HI]. First three-dimensional (3D) simulations 
based on this model were published in [15j, showing 3D 
resistive filamentation. Grcmillct ct al. also derived the 
linear growth rates, now depending on resistivity rj. The 
present version of the model is described in more detail 
in 0. 

Recently, two-dimensional simulations of the cone- 
guided target experiment [a, Ml were published by Camp- 
bell et al. [17J and Mason [l8j . Based on different hybrid 
codes, they could reproduce the measured core heating 
of 800 eV. Using the present model, we have obtained 
similar 2D results comparisons with the linear the- 
ory of the resistive filamentation instability are found in 
[20j. Here we present first three-dimensional simulations 
of electron transport and deposition in the high-density 
part of a fast ignition target. It is shown how 3D beam fil- 
amentation seeded in the corona leads to multi-hot-spot 
heating of the core close to DT ignition. The present 
study is motivated by the next generation of high-power 
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FIG. 1: Central cut through imploded target configuration: 
(a) isocontours of density in g/cm 3 , (b) density profile at y=75 
fim. 



facilities to demonstrate fast ignition 21, 2^ . 

Simulation parameters. The simulated configuration is 
shown in figure Q] It consists of 0.2 mg DT fuel com- 
pressed into a spherical blob of 400 g/cm 3 peak den- 
sity and 75 jum diameter (full width at half maximum, 
FWHM); it sits on a density pedestal of 1.5 g/cm 3 (the 
halo). A beam of fast electrons is injected from the left at 
2 = 0. We imagine that it emerges from the tip of a cone 
at this position. The cone itself and the laser pulse gen- 
erating the beam inside the cone are not simulated here. 
We rather model the injected beam in form of directed 
Gaussian electron distributions in radius and time with 
a spot radius of 20 fim and a duration of 10 ps, both at 
FWHM. The pulse has a power of 6 PW, a total energy 
of 60 kJ («30 kJ within FWHM) and is centered at 7 ps. 
The energy distribution of the beam electrons is assumed 
to be TD relativistic Maxwellian with temperatures de- 
pending on the local laser intensity /, also assumed Gaus- 
sian in radius and time, by the ponderomotive scaling 
formula T b » / • m e c 2 [(1 + JA 2 /13.7 GW) 1 / 2 - 1], where 
A is the wavelength. PIC simulations 23j give front fac- 
tors / ~ 1 — 3, depending on the scale-length of the 
plasma in which the electrons are accelerated. For cone- 
guided fast ignition with 10 ps pulse durations and elec- 
tron acceleration along the cone surface 2J] , the factor is 
expected to be larger than / ~ 1 which applies to sharp 
surfaces. Here we consider different cases with mean elec- 
tron kinetic energies (averaged over the FWHM of the 
distributions in radius and time) in the range of (E) — 
1.5 - 5.5 MeV; 2.5 MeV is taken as a reference value. 
This mean energy corresponds to a laser irradiance of 
1.5 x 10 20 W/cm 2 (FWHM) at 0.35 fim, assuming laser- 
to-electron transfer of 50% and beam compression by a 
factor of 3 due to geometrical cone convergence, in agree- 
ment with the results reported in 0, 0, Hj, [25| . The ini- 
tial angular distribution of fast electrons is obtained as 
in (2d| . Electrons with energy E = (7 — l)m e c 2 are in- 
jected with a randomly chosen half-angle between and 
tan _1 [/i-\/2/(7 — 1)]. The parameter h is used to adjust 
the initial beam opening half-angle as 22.5° (FWHM), 
consistent with the cone experiment @, 0] and the simu- 
lations in [13, El- 
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FIG. 2: (a) Central longitudinal cut of magnetic field B x in 
kT generated by electrons with a mean kinetic energy of 2.5 
MeV at the peak of the pulse, (b) perpendicular cut of beam 
current density j z in units of 10 14 A/cm 2 at z=98 fim. Plasma 
densities higher than 200 g/cm 3 are located inside the dashed 
circle. 



The imploded fuel configuration shown in figure [I] has 
been scaled from that reported in 



17, m 



The 



parameters of this configuration, the peak density of 400 
g/cm 3 , the distance of 150 fim from cone tip, and the ini- 
tial plasma temperature, have been adapted from implo- 
sion simulations of cone targets with the code SARA-2D 
27J ■ The plasma resistivity depends on the temperature 
distribution. The SARA-2D simulations indicate temper- 
atures in the range of 300 eV to 1 keV. For simplicity, a 
uniform initial DT temperature of 500 eV is taken here, 
which sets the initial resistivity to a level of 10~ 8 fim. 
Concerning the numerical parameters, we have chosen a 
cell width of 1 fim in each coordinate, a time step of 
3 fs, and a total number of 3.6 x 10 7 particles injected 
over the time interval of - 14 ps. Free boundaries have 
been used in all simulations. Classical Spitzer resistivity 
is chosen for the DT plasma, and MPQeos tables [28| are 
used to compute electron and ion temperatures from the 
deposited energy. 

Results. The injected current of 3.5 GA decays into 
filaments after a propagation distance of z = 70 ^m; 
this is seen in figure [UJa) in terms of the B-field and 
in figure [2](b) in terms of current density. Actually we 
find that the filaments start to grow in the halo region 
(z < 50 fim) and are then strongly amplified in the den- 
sity slopes of the blob. The growth rate is consistent with 
the linear theory developed in 15j. Resistive filamenta- 



tion scales with plasma resistivity and therefore depends 
strongly on electron temperature T e . In the lower density 
region, ohmic heating and Coulomb energy deposition 
lead to high electron temperatures with a mean value of 
(T e ) rs 40keV, much higher than the ion temperature Tj, 
and here magnetic fields saturate at levels of 100 T due to 
low resistivity. At higher densities (z > 70 fim), sufficient 
energy transfer from electrons to ions takes place, and we 
find T e ss Ti < 20keV. The corresponding ion tempera- 
ture is plotted in figure[3](a). Here higher resistivity leads 
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FIG. 3: Ion temperature of DT in units of log 10 (Ti/eV) at the 
end of the pulse of electrons with a mean kinetic energy of 2.5 
MeV. (a) longitudinal cut at x = 75 /im, (b) transverse cut 
at z = 120 fim. Beam-generated fields have been artificially 
suppressed in (c) and (d). Plasma densities higher than 200 
g/cm 3 are located inside the dashed circle. 



to .B-field of rj 1 kT and enhanced filament growth. One 
should notice that the resistive filamcntation observed 
here is weaker than that obtained from collisionless PIC 
simulations 
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12| and has a wider spatial scale. Given 
the broad energy spectrum of beam electrons only those 
with lower energies (< 1.5 MeV) are trapped in the cur- 
rent channels, while the others scatter freely and tend to 
smooth the structure. We find that both temperatures 
and £?-fields are well described by the analytic scaling 
laws derived by Bell and Kingham k 2j|. We have also 
checked that filamentation persists for initial tempera- 
tures of 1 keV and halo densities up to 10 g/cm 3 , but 
disappears for beam kinetic energies (E) > 4.5 MeV. 

The filaments shown in figure [2jb) carry about 10 MA 
beam current each. This beam current is almost com- 
pletely compensated by the plasma return current, im- 
plying j r rj —jb', the net current of a filament is only 
about 10 kA, consistent with the magnetic field strengths 
of 1 kT observed in figure [Ifa) . It is worthwhile em- 
phasizing that the filamented current distribution heats 
electrons significantly more than a uniform current dis- 
tribution of same total current because of the j% depen- 
dence of ohmic heating. This is how filamcntation leads 
to enhanced beam stopping within the present model. 
It differs from the strong anomalous stopping found in 
PIC simulations [HI, [HI for lower plasma densities (about 
10 x beam density). Sentoku et al. have interpreted this 
stopping as stochastic scattering of return current elec- 
trons by magnetic perturbations \B\ giving rise to an ad- 
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FIG. 4: Pressure of DT in Tbar for the electron pulse with 
a mean kinetic energy of 2.5 MeV: (a) longitudinal cut at 
x = 75 fim and (b) transverse cut at z = 120 /im. Plasma 



densities higher than 200 ; 
circle. 



/cm are located inside the dashed 



ditional effective resistivity ?7e// ~ \B\/(en p ). Even tak- 
ing the scaling |B| oc n^ 5 from j2j|, rj e ff decreases with 
plasma density n p . For the parameters of the present 
simulation, we find f] e ff to be always smaller than the 
Spitzer resistivity. Therefore, we conclude that anoma- 
lous stopping of this kind plays no significant role here, 
in agreement with the results of Mason [l8| ■ 

Figures [3] and [4] present the central results of this pa- 
per, showing DT fuel close to ignition. Notice that beam 
filamentation in the corona is responsible for the frag- 
mented heating pattern in the core. Artificial suppression 
of the beam-generated fields (i.e. Coulomb deposition 
only) would lead to smooth core heating with a maximum 
temperature of 4 keV. This is shown in Figs. [3jc) and 
(d), for comparison. With fields present, a multi- hot-spot 
ignition region is formed in the high-density fuel with 
maximum temperatures beyond 10 keV. The conjecture 
here is that this will help ignition due to the nonlinear 
scaling of fusion reactivities with temperature. Actually, 
this needs to be confirmed in more detailed simulations, 
including hydrodynamics and fusion reaction physics in 
3D geometry. 

Here we give some estimates, based on the pressure 
distribution shown in figure [fj The DT ignition condi- 
tion (see page 85 in can be written in compact form 
asphRh > 45{Phl Pc) 1 ^ 2 Tbar/mi, where index h refers to 
hot fuel and index c to surrounding cold fuel. This condi- 
tion combines the threshold values for phRh and Th and 
holds for temperatures 5 < T^/keV< 15. From Figs. 02a) 
and (b), we find for the central hot spot 2phRh ~ 50 
Tbar ^m in longitudinal direction. In transverse direc- 
tion, it is 2phRh ~ 15 Tbar /xm, but here a number of 
neighboring hot spots will cooperate and 2p c R c w 50 
Tbar /xm, obtained from p c = 1 Tbar and 2R C = 50 /im, 
may serve as an estimate. We conclude that the refer- 
ence case shown in Figs. [3] and |4] is close to ignition. It 
should be understood that the core heating is almost ex- 
clusively due to Coulomb deposition of beam electrons. 
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Ohmic heating by return currents dominates in the halo, 
but plays only a minor role for the overall energy bal- 
ance. Beam-generated fields turn out to contribute to 
the core heating indirectly, mediated by filamentation 
and collimation effects. Beam collimation is observed 
in figure [3fa) when compared with figure G3c). 

Figure [5] shows what fractions of the injected beam 
energy are deposited in different parts of the target. The 
total deposition, given versus beam kinetic energy (E), 
drops from 90% to 40%, when raising (E) from 1.5 to 5.5 
MeV. Most of the deposition is due to classical Coulomb 
collisions, consistent with the areal density of 2.9 g/cm 2 
along the axis. Notice that the other part of the energy is 
not deposited at all, but passes through the target and is 
lost. Clearly this makes average beam energies beyond 5 
MeV prohibitive. The very important partition between 
deposition into high density core and low-density zones 
with p < 200 g/cm 3 is also shown in figure \5[ We find 
that the energy coupling to the core amounts to 30% at 
1.5 MeV and 20% at 5 MeV. It is less sensitive to (E) than 
the coupling into the lower density regions. Of course, the 
core coupling strongly depends on the divergence angle 
of the injected beam, which is therefore a key parameter 
for fast ignition. In the present simulations, the core 
coupling efficiency degrades by 40%, when raising the 
angle from 22.5° to 30°. On the other hand, one should 
notice that magnetic pinching of the beam improves core 
coupling significantly. Suppressing all beam-generated 
fields in the simulations would lead to the dashed curves 
in figure El 

Estimating the laser-to-fast-electron conversion effi- 
ciency to be 50%, we find a laser pulse energy of 100 
- 150 k J necessary to ignite a target. There may be pos- 
sibilities to reduce this energy, e.g. by shortening the 
distance between cone tip and blob or b y ca reful design 
of the cone to reduce beam divergence [30(. Certainly, 
transport in the cone needs to be included in more com- 
plete studies, in particular to account for the potential 
barrier [l8[ and filamentation 24 1 at the cone tip. 

In conclusion, the message of this paper concerning fast 
ignition of inertial fusion targets is that a giga-ampere, 
multi-PW current can be transported through the steep 
gradients of the plasma corona toward the high-density 
fuel core. This is shown here for the first time in 3D 
geometry, using hybrid PIC simulation. Central ques- 
tions could be answered: Collective magnetic effects play 
a major role for core heating, but in an indirect way. Re- 
sistive beam filamentation grows in the low-density halo 
and seeds the 3D multi-prong beam, which then pene- 
trates the core. Of course, 3D simulation is crucial in 
this context. 

In the core, collective behavior is suppressed due the 
large plasma-to-beam density ratio, and energy deposi- 
tion takes place almost exclusively by classical Coulomb 
collisions. We find a fragmented hot spot configuration, 
and the fragmentation may actually help fuel ignition, 
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FIG. 5: Fraction of total pulse energy deposited in the target 
(squares), in low-density (<200 g/cm 3 ) zones (triangles) and 
in high-density (>200 g/cm 3 ) core (diamonds). Solid lines 
correspond to full simulations with beam-generated fields and 
Coulomb energy deposition. Dashed lines correspond to sim- 
ulations with beam-generated fields suppressed. The fraction 
due to ohmic heating is small throughout. 



since concentrating the energy in a number of prongs 
rather than heating the whole volume spanned by the 
prongs leads to higher temperatures and therefore to 
more heating by fusion products. The (pR) values ob- 
tained for the reference case indeed indicate that this case 
is close to ignition. More detailed simulations including 
hydrodynamics and fusion heating are now in progress to 
confirm this point. 

Concerning collective beam deposition, we find, within 
the physical model used here, that indeed beam fila- 
mentation enhances ohmic heating, because it depends 
quadratically on the return current density j r . But this 
additional deposition is not identical with the anoma- 
lous stopping found in PIC simulations at lower densi- 
ties. It also contributes little to the overall energy bal- 
ance in the present simulation of fast ignition. Rather 
the self-generated B-fields help by collimating the rela- 
tivistic beam, and this improves the coupling efficiency 
substantially. 
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